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RESUMEN 

El resumen sera traducido al espanol por los editores. We present the first results obtained from 
our pilot ultra-deep near-IR survey with ISAAC/VLT, aimed at the detection of z> 7 sources using lensing 
clusters as natural Gravitational Telescopes. Evolutionary synthesis models of PopIII and extremely metal-poor 
starbursts have been used to derive observational properties expected for these high-z galaxies, such as expected 
magnitudes and colors, line fluxes for the main emission lines, etc. These models have allowed us to define fairly 
robust selection criteria to find z^ 7—10 galaxies based on broad-band near-IR photometry in combination 
with the traditional Lyman drop-out technique. The magnification in the core of lensing clusters improves the 
search efficiency and subsequent spectroscopic follow up. The research efficiency will be significantly improved 
by the future near-IR multi-object facilities such as EMIR/GTC and KMOS/VLT-2, and later the JWST. 

ABSTRACT 

We present the first results obtained from our pilot ultra-deep near-IR survey with ISAAC/VLT, aimed at 
the detection of z> 7 sources using lensing clusters as natural Gravitational Telescopes. Evolutionary synthe- 
sis models of PopIII and extremely metal-poor starbursts have been used to derive observational properties 
expected for these high-z galaxies, such as expected magnitudes and colors, line fluxes for the main emission 
lines, etc. These models have allowed us to define fairly robust selection criteria to find z^ 7—10 galaxies 
based on broad-band near-IR photometry in combination with the traditional Lyman drop-out technique. The 
magnification in the core of lensing clusters improves the search efficiency and subsequent spectroscopic follow 
up. The research efficiency will be significantly improved by the future ncar-IR multi-object facilities such as 
EMIR/GTC and KMOS/VLT-2, and later the JWST. 
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1. MOTIVATION 

The detection of the very first stars and galaxies 
forming from pristine matter in the early Universe 
remains one of the major challenges of observational 
cosmology (cf. review of Loeb & Barkana 2001; Weiss 
et al. 2000; Umemura & Susa 2001). These sources 
constitute the building blocks of galaxies, forming 
at redshifts z ~ 6-30. Indeed, the recent WMAP 
results seem to place the very first generation of 
stars at z ~ 10-30 (Kogut et al. 2003). Important 
progress during the last years has permitted direct 
observations of galaxies and quasars out to redshifts 
z ~ 6.6 (Hu ct al. 2002, Fan et al. 2003, Kodaira et 
al. 2003). At z> 7, the most relevant signatures are 
expected in the near-IR A> 1/im window. Modeling 
efforts during the last years have been motivated by 
JWST, which should be able to observe these ob- 
jects at redshifts up to z ~ 30. However, the delay 
of JWST and the availability of well suited near-IR 
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facilities in ground-based 8-10m class telescopes has- 
ten the development of observational projects target- 
ting z > 7 sources. The first studies on the physical 
properties of these sources can be presently started 
with instruments such as ISAAC/VLT, and contin- 
ued with the future multi-object spectrographs such 
as EMIR at GTC (~ 2006), or KMOS for the second 
generation of VLT instruments (~ 2009). 

In this paper we present the first results from 
our ongoing pilot project with ISAAC/VLT, aiming 
at looking for the first stars/galaxies through lens- 
ing clusters used as Gravitational Telescopes (GTs). 
The plan of the paper is the following. We sum- 
marize in section 2 the observational properties ex- 
pected for galaxies at 7 < z < 10. In section 3 we 
present the main characteristics of our ISAAC /VLT 
project. The method used to identify sources and 
the first results obtained are also presented in this 
section, in particular the first spectroscopic confir- 
mation of a z=10 candidate. The implications of 
these results for EMIR/GTC are briefly discussed in 
section 4. 
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2. PROPERTIES OF 7 < Z < 10 GALAXIES 

We have used the evolutionary synthesis models 
by Schaerer (2002, 2003) for Population III and ex- 
tremely metal deficient stars to derive the expected 
magnitudes and colors of galaxies at 7 < z < 10. 
The dependence of their properties on the IMF and 
upper mass limits for star formation were studied in 
previous papers (e.g. Schaerer & Pello 2001; Pello & 
Schaerer 2002). We have also computed the expected 
S/N ratios for the main emission lines using tele- 
scope and instrumental parameters in the near-IR 
corresponding to ISAAC/ VLT, and the future multi- 
object facilities EMIR/GTC and KMOS/VLT. 

According to our simulations, the predicted mag- 
nitude in the K'(Vega) band is typically ~ 24 to 25 
for the reference stellar halo mass of 10 7 M . Neb- 
ular continuous emission and strong emission lines 
have important effects on the integrated fluxes and 
colors of such galaxies. Figure 1 displays the J-H 
versus H-Ks color-color diagram for different mod- 
els with 5 < z < 11, compared to the location of 
stars and normal galaxies at < z < 8, including 
the starbursts templates SB1 and SB2, from Kinney 
et al. 1996, and the low metallicity galaxy SBS0335- 
052. Broad-band colors do not allow to constrain 
the physical properties of these galaxies, such as the 
IMF, star-formation mass range, ages, etc., but they 
are useful to identify sources on ultra-deep photo- 
metric surveys (cf. below). 

The main emission lines expected are Lyman a 
and Hell A 1640, the strongest Hell line. In princi- 
ple, Lyman a can be detected on near-IR spectro- 
scopic surveys with 8-10 m class telescopes, with a 
reasonable S/N over the redshift intervals z ~ 7 to 18, 
with some gaps depending on the spectral resolution 
(OH subtraction), the atmospheric transmission, the 
intrinsic properties of sources and the intergalactic 
medium (IGM) transmission. A joint detection with 
Hell A1640 should be possible for z ~ 5.5-7.0 (Ly- 
man a in the optical domain), z ~ 7-14 with both 
lines in the near-IR. A measure of the hardness of 
the ionising flux, which constrains the upper end of 
the IMF and the age of the system, could be ob- 
tained through the detection of both Hell A1640 and 
Lyman a. A rough measurement of the continuum 
using the photometric spectral energy distribution 
(SED) provides additional information on the stellar 
populations, extinction, etc. 

3. A PILOT PROJECT WITH ISAAC/VLT 

Lensing clusters acting as GTs are useful tools to 
investigate the properties of distant galaxies. Suc- 
cessful systematic searches for z > 5 galaxies with 



All models 



X 



A \# ! 






\ ¥^>-' 










\ f 

- JSf i ? 


. • . bf \m — — 

JbiSbSX Ijiy^ 'r " 0<z<8 normal galaxies 










starS" "I "^t* 


f Kinney Bursts ^7 
| SBS0335 / 

j Reddening 





-0.5 0.5 1 1.5 3 2.5 3 3.5 



H-K 

Fig. 1. J-H versus H-Ks color-color diagram (Vega sys- 
tem) for young starbursts (metal-poor and metal-free) 
objects, showing the position for the different redshifts 
over the interval of z ~ 5 to 12. The position of stars 
and normal galaxies up to z < 8 including also starbursts 
templates (SB1 and SB2, from Kinney et al. 1996, and 
the low metallicity galaxy SBS0335-052) are also shown, 
as well as the shift direction induced by reddening. Var- 
ious line styles (solid, dashed, dotted) illustrate young 
starburst models with Lyman-a contributions of 100, 50, 
and %. Because of the limiting magnitudes in the differ- 
ent filters, most z ~ 7-10 candidates are expected within 
the "High-z candidates" region. This diagram illustrates 
that such galaxies can well be separated from "normal" 
objects on the basis of JHK photometry. 

GTs have appeared during the last three years (e.g. 
Ellis et al. 2001; Hu et al. 2002; Kneib et al. 2004). 

We were granted ISAAC/VLT time to develop 
a project aimed at searching for z > 7 galaxies us- 
ing GTs. This prototype observational program is 
presently going on. The method used and the first 
results obtained are summarized below. 

3.1. Method 

Galaxy candidates at z > 7 are selected from 
ultra-deep JHKs images in the core of gravitational 
lensing clusters for which deep optical imaging is also 
available, including HST data. We apply the Lyman 
drop-out technique to deep optical images. The main 
spectral discontinuity at z > 6 shortward of Lyman 
a is due to the large neutral H column density in 
the IGM. In addition, galaxies are selected accord- 
ing to Fig. 1: we require a fairly red J — H color (due 
to the discontinuity trough the J band), and a blue 
H — Ks colour corresponding to an intrinsically blue 
UV rcstframe spectrum. The detection in at least 
two bands longward of Lyman a and the combina- 
tion with the above H — Ks colour criterion allows 
us to avoid contamination by cool stars. 
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Fig. 2. The object #1916 in A1835: a) 2D spectrum 
showing the detected emission line at 1.33745 /im(leading 
to z — 10.00 when identified as Lyman-a, as well the 
nearby field galaxy and the [O ill] A5007 line of the galaxy 
#2582 (z — 1.68, Richard et al. 2003) used as a reference 
to stack the data sets. The line is seen on 2 independent 
overlapping bands, b) Photometric optical to near-IR 
SED in Fx units, c) and d) Thumbnail images in the 
HST-WFPC2 R band, where the object is not detected, 
and in the H band. 

Ultradeep JHKs imaging of 2 lcnsing clusters 
were obtained: A1835 (z = 0.25) and AC114 (z = 
0.31), under excellent seeing conditions (~ 0.4" to 
0.6"). A complete report and analysis of these ob- 
servations will be given elsewhere (Richard et al. 
in preparation). The best candidates are retained 
for follow-up spectroscopy with ISAAC/VLT. The 
broad band SED of each source is used to constrain 
the redshift using a modified version of the public 
photometric redshift code Hyperz (Bolzonella et al. 
2000), including numerous empirical and theoretical 
spectral templates for both galaxies and AGNs, and 
keeping internal extinction as a free parameter. The 
J band domain, where Lyman-a should be located 
for objects in this redshift interval, is systematically 
explored with ISAAC according to the priorities set 
by the photometric redshift probability distribution. 

3.1.1. First Results on 7 < z < 10 Photometric 
Candidates 

Applying the above broad-band search criteria 
has provided 5-10 z > 7 galaxy candidates per lens- 



ing cluster in the observed fields. About half of 
them have photometric SEDs accurate enough to al- 
low a robust identification within this redshift in- 
terval, and thus to be included in the spectroscopic 
follow up sample. Lensing effects are taken into ac- 
count to determine the intrinsic properties of these 
sources and the effective areas surveyed on the dif- 
ferent source planes. Typical magnification factors 
range between > 2 and > 10. The corresponding 
number density of objects within 7 < z < 10 ranges 
between 2 x 10~ 2 and 4x 10~ 4 per Mpc -3 , i.e. typ- 
ically a few 10 3 objects deg -2 . The main uncertain- 
ties are due to simplifying assumptions, incomplete- 
ness corrections, and the values adopted for the typi- 
cal magnification factors of our sample. From the ob- 
served magnitudes and magnification factors derived 
from the lens model their star formation rates are es- 
timated to be a few M Q yr _1 (lcnsing corrected, as- 
suming a "standard" Salpeter IMF from 1-100 M Q ). 
Adopting typically SFR ~ 2 MQyr" 1 we obtain an 
estimate of the star formation rate density shown 
in Fig. 3, and compared to previous estimates be- 
tween z = — 7. Obviously for now the SFR density 
at z > 7 remains highly uncertain. Ricotti et al. 
(2004) recently derived a fairly high SFR density at 
z = 10 based on our galaxy detection with spec- 
troscopic redshift (A1835 IR1916, cf. below). How- 
ever, their study suffers mostly from a highly uncer- 
tain survey volume area and an overestimated mean 
magnification factor. In short, a more detailed anal- 
ysis and additional observations are required to draw 
more firm conclusions on the behaviour of the SFR 
density at z > 7. 

3.1.2. Spectroscopic confirmation of a z=10 
candidate 

We report in a recent paper the first likely spec- 
troscopic confirmation of a z = 10.0 galaxy in our 
sample (see Pello et al. 2004). This galaxy (called 
#1916) was the best candidate found in the field 
of the lensing cluster A1835, and a good example 
of the research procedure. In this case, the photo- 
metric redshift probability distribution shows a clear 
maximum at redshift z p hot ~ 9-11. This estimate is 
mainly corroborated by a strong break of > 3.1 to 
3.7 AB mags between VRI and H, which has a high 
significance independently of the definition used for 
the limiting magnitudes (see Fig. 2). 

The spectroscopic redshift determination was ob- 
tained with ISAAC/VLT between 29 June and 3 July 
2003, under excellent seeing conditions, using a 1 arc- 
sec slit width. The observed spectral interval covers 
the range from 1.162 to 1.399 /mi(redshifts z ~ 8.5- 
10.5 for Lyman-a:). The observations resulted in the 
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Fig. 3. SFR density derived from our 5 best candidates 
at 7 < z < 10 in the lensing cluster A1835, compared to 
previous estimates up to z < 7. Upper and lower limits 
correspond to extreme assumptions for the typical mag- 
nification factors ranging between > 2 and > 10. 



detection of one weak emission line at the 4-5 a level 
with an integrated flux of (4.1 ± 0.5) x 10~ 18 erg 
cm~ 2 s _1 at a wavelength of 1.33745 fan (see Fig. 2), 
which appears on 2 different overlapping wavelength 
settings. When identified as Lyman-a in good agree- 
ment with the photometric SED, the observed line 
gives a redshift of z = 10.00175 ± 0.00045, the most 
likely one given the data set (see a detailed discus- 
sion in Pello et 2004). 

Given the location of the image on the lensing 
plane, the source is highly magnified by a factor of 
25. Its intrinsic AB magnitude is 28.5 and 29 in 
the H and Ks bands respectively. The star for- 
mation rate obtained from a direct scaling of its 
UV restframe flux ranges between 1.8 and 3 M Q 
yr -1 after lensing correction, assuming a "standard" 
Salpeter IMF from 1-100 M Q . The UV restframe 
spectrum is very blue, compatible with a small ex- 
tinction in this galaxy. Only a relatively small frac- 
tion of the Lyman-cv flux emitted from the source is 
observed. Scaling young burst models with metal- 
licities Z > 1/50Z Q to the observed SED we obtain 
A/* ~ 2 x 10 s M (or lensing corrected values of 
A/* - 8 x 10 6 M ) and luminosities of L ~ 4 x 10 11 
Lq (2 x 10 10 Lq lensing corrected) for the standard 
Salpeter IMF. The observed and derived properties 
of this object are in good agreement with expecta- 
tions for a young protogalaxy experiencing a burst 
of star formation at z = 10. 



4. IMPLICATIONS FOR EMIR/GTC 

Simulations and pilot observations conducted by 
our group on the possible detection of z > 7 galax- 
ies support this scientific case in view of the fu- 
ture ground-based near-IR facilities such as EMIR at 
GTC and the future KMOS for the VLT 2nd genera- 
tion (Schaerer & Pello 2001; Pello & Schaerer 2002). 
An important issue of this project in view of the fu- 
ture facilities is the strategy for target selection, and 
the efficiency achieved in spectroscopic studies. As 
shown in sections 2 and 3.1, starbursts could be 
detected from deep near-IR photometry based on a 
measurement of two colors with accuracies of the or- 
der of 0.3-0.5 mag. As shown in Fig. 4, GTs are 
the ideal fields for the first prospective studies. Wc 
compare in this figure the expected number counts of 
high-z galaxies with K'< 24 in our GT ongoing sur- 
vey, for 2 different extreme assumptions for the IMF, 
leading to intrinsically bright or faint sources. A sim- 
ple Press-Schechter formalism for the abundance of 
halos and standard ACDM cosmology were consid- 
ered in these order of magnitude simulations, as well 
as a conservative fixed fraction of the baryonic mass 
in halos converted into stars. 

As seen in Fig. 4, strong lensing fields are a fac- 
tor of ~ 5 — 10 more efficient than blank fields of the 
same size in the z ~ 8 — 10 domain. This is the moti- 
vation of all present searches for z > 6-7 galaxies using 
GTs (and ISAAC in particular). The large field of 
view of the future instruments (between ~ 6' — 8' for 
EMIR/GTC and KMOS/VLT) provides the same ef- 
ficiency than GT fields in "optimistic" models, and 
even better efficiency up to z~10, simply because the 
size of the field compensates the magnification (and 
dilution) effects. In the most pessimistic models (as- 
suming intrinsically fainter sources), the benefit of 
the wide field of view versus strong lensing extends 
only up to z~7.5, the magnification being a major 
benefit in this case. As shown in the figure, when 
this regime is attained, the number of sources ex- 
pected per redshift bin Az = 2 falls below 1 source 
per arcmin , and thus we can still take benefit from 
the presence of GTs in the surveyed field. The red- 
shift bin considered (Az = 2) corresponds to the 
typical wavelength interval which can be explored in 
a single shot with EMIR. These numbers are to be 
considered as orders of magnitude. Using higher red- 
shift GTs could improve the situation for " compact" 
lensing studies with respect to ~ 6' — 8' blank fields. 

Thanks to its multiplexing capabilities, spectral 
resolution and wide field of view, EMIR at GTC will 
be the best suited instrument to continue exploring 
the formation epoch of the first stars in galaxies us- 
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Fig. 4. Comparison between the expected number counts 
of "primordial galaxies" up to K<24, per redshift Az = 
0.1, with and without the presence of a strong lensing 
cluster, AC114 (z=0.31) in this example, within the field 
of view of ISAAC/VLT, compared to EMIR/GTC. Two 
different extreme assumptions for the IMF are used. The 
solid horizontal line displays the domain corresponding 
to less than 1 source per arcmin -2 within Az = 2. 

ing significant samples of high— z sources. A rela- 
tively high spectral resolution (R ~ 3000-5000) con- 
siderably increases the chances of detection between 
the sky lines. Once this is obtained, the data can 
be rebinned to increase the S/N. The medium spec- 
tral resolution is also favoured to attempt to measure 
the emission line profiles, in particular to distinguish 
genuine "primordial" stellar sources from potential 
very high-z AGN (cf. Tumlinson et al. 2001). 
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